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TRANSONIC  SHOCK  INTERACTION  WITH  A 
TANGENT I ALLY- INJECTED  TURBULENT  BOUNDARY  LAYER 

4.  4»  4- 

G.  R.  Inger  and  A.  Deane 
West  Virginia  University,  Morgantown,  W.  Vj. 


Abstract 


A  non- asymptotic  triple  dec<  theory  oi  tran¬ 
sonic  shock  turbulent  boundary  layer  interaction  ** 

is  described  which  takes  into  iccoant  the  influence  *  . 
of  a;  stream  tangential  injection  on  a  curved  wall. 

In  iddition  to  rvjm.rer  and  the  cl 

T 

strength,  the  thoorv  is  p.ir  tr.otoriroa  by  iroitrarv  c  1 

T 

values  >f  the  inco^m;  bound  try  layer  snipe  t'-etor, 
wall  jet  maximum.  velocity  rati  and  the  non- 
1  i" ens ionu 1  height  of  this  rati  ';  results  of  a 
comprehensive  pur  i-'etr  ic  study  v»re  then  ['resented. 

It  ic  shown  that  the  wall  jet  effects  significantly 
reduce  i'oth  the  -treamwise  scale  and  d  ispi  ent 
thickening  of  the  interaction  zone.  While 
increasing  the  upstream  and  downstream  skip,  fric- 
tion  levels,  these  effects  also  rv  iuce  the  p  in imam 
interactive  Of  and  thus  hasten  the  onset  of 
incipient  separation  at  the  shock  foot. 


specific  heat  ratio 
boundary  layer  thickness 
boundary  layer  displacement  thickness 
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Nomenclature 

Yin  Dr iest-Ce‘>eci  wall  turbulence  :  imping 
I'arameter 

skin  friction  coefficient,  2 x 

e  e 
o  o 

skin  f  r  ic t  ion  i nc r om.ent  due  t o  wa  11  jet 

pressure  coefficient,  2  p'  .  U 

e  0 
o  o 

1  it^*r  «l  spreading  factor  of  wall  i^t 
boundary  lay^r  sh  ipe  factor,  5*  * 

incompressible  shape  factor 
curvature  of  wall  in  interaction  region 
Mach  number 
static  pressure 

interactive  pressure  per  turb.it  ion  ,  p-p^ 
pressure  jump  across  incident  shock 
Reynolds  numbers  biuod  on  length  .  and 
hound  u  y  layer  thickness,  respectively 
non-dimensional  wall  shear  function  of 
wall  jet 

absolute  temperature 
h  is ic  interactive  wail-turbulence 
par ameter 

s^reamwise  and  normal  interactive  distur¬ 
bance  velocity  components ,  respectively 
wall  jet  component  of  total  velocity 
;  -rot  i  le 

undisturbed  incoming  boundary  layer 
velocity  in  x-direction 
streamwise  and  normal  distance  coordin¬ 
ates  (origin  at  the  inviscid  shock  in¬ 
tersection  with  the  wall) 
effective  wall  shift  seen  by  interactive 
incise id  flow 
Location  of  YU 

max 


•’rofensor  »nd  Associate  Chairman,  Dept,  of 
^Mechanical  and  Aerospace  Engineering 
Graduate  student. 


Subscrip t  s 


AD  adiabatic  wall 

1  undisturbed  inviscid  values  ahead  of 

incident  shock 

e  conditions  at  the  boundary  layer  edge 

me  incompressible  value 

inv  inviscid  disturbance  solution  value 

max  velocitv  profile  maximum  due  to  wall  jet 

a  undisturbed  incoming  boundary  layer  pro¬ 

perties 


1 .  Introduct ion 


The  use  of  tangential  slot- in ject ion  to  influ¬ 
ence  and  control  turbulent  boundary  layer  behavior 
has  been  extensively  studied  in  various  types  of 
lowspeed  external  and  internal  aerodynamic  flow 
fields  (e.  g.f  on  circulation-controlled  airfoils, 
slotted  slaps,  in  film  cooling  applications  and 
for  separation  control  in  inlets  and  diffusers). 

»n  recent  years,  many  applications  of  such  injec¬ 
tion  have  arisen  in  supercritical  transonic  flow 
fields  where  local  shock  wave  is  present;  however, 
little  is  presently  available  to  provide  a  basic 
understanding  of  how  the  resulting  shock-boundary 
interaction  (’'SBLI")  alters  the  influence  of 
tangential  injection.  Conversely,  in  such  super¬ 
critical  flows  it  may  be  of  interest  to  know  how 
the  effects  of  SBLI  may  be  altered  by  the  use  of 
injection.  The  present  paper  addresses  these 
questions  for  the  case  of  steady  non-separating  2-D 
turbulent  boundary  layers  on  adiabatic  surfaces  of 
small-to-moderate  longitudinal  curvature. 

The  primary  objectives  of  our  work  are  to 
develop  a  fundamental  theory  of  a  transonic  SBLI 
region  occurring  downstream  of  a  tangentially- 
injected  turbulent  boundary  layer  on  a  curved  wall 
{Fig.  1)  and  then  to  present  the  results  of  a 
parametric  study  of  this  theory  showing  the 
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relat ionship  between  the  dominant  physical  parame¬ 
ters,  the  injection  and  the  physics  of  the  SBLI 
zone.  In  Section  2,  we  briefly  outline  the  non- 
asymptotic  triple  deck  theory  of  a  S3LI  zone  on  a 
curved  surface  without  tangential  injection.  Then 
by  taking  the  SBLI  zone  sufficiently  far  down¬ 
stream  of  the  injection  slot  for  mix  in}  of  the 
wall  jet  and  overlying  turbulent  boundary  layer  to 
have  produced  a  well-defined  "jet-bulged"  boundary 
layer  profile,  the  interactive  perturba¬ 

tion  fiel'1  c  tused  by  normal  shock  interaction  with 
this  profile  1=  analyzed  in  Section  3  by  an  exten¬ 
sion  of  the  aforementioned  SBLI  theory.  This  is 
followed  in  Section  4  by  presentation  and  discus¬ 
sion  of  the  results  of  a  parametric  study  of  this 
extended  solution  for  the  interactive  pressure, 
displacement  thickness  and  skin  friction  effects. 


profile  model  due  to  Walz^,  which  is  characterized 
not  only  by  the  shock  Mach  numbe.  Mj  and  the  boun¬ 
dary  layer  thickness  Reynolds  number  Ro^#  but  -also 
by  arbitrary  nonequilibrium  values  of  the  incom¬ 
pressible  shape  factor  The  resulting  predic¬ 

tions,  such  as  typically  illustrated  in  Fig.  3, 
show  that  H i ^  has  a  very  large  effect  on  the  local 
and  downstream  interactive  properties  that  is 
important  to  account  for  in  practical  applications 
By  thereby  accomodating  a  wide  range  of  possible 
upstream  histories  of  pressure  gradient,  heat  and 
mass  transfer,  the  theory  has  found  wide-spread 
success  as  an  interactive  m.odule  in  global  compos¬ 
ite  viscous-inviscid  flow  field  analysis  programs 
on  supercritical  airfoils  and  projectiles,  while 
also  proving  adaptable  to  the  accomodation  of  new 
ef  f  ects . 


2.  Brief  Outline  of  the  Basic  SBLI  Theory 


2 . 1 1  The  Tr t pie -Deck  Model 

It  is  well-known  experimentally  that  when 
separation  occurs,  the  disturbance  flow  pattern 
associ  with  normal  shock-boundary  layer  inter¬ 
action  is  a  very  complicated  one  involving  a 
bifurcated  shock  oattern*,  whereas  the  un.^-pa rated 
case  pertaining  to  turbulent  boundary  layers  up  to 

1.3  has  instead  a  much  simpler  type  of  inter¬ 
action  pattern  which  is  more  amenable  to  analytical 
treatment  (Fig.  2).  The  flow  consists  of  a  known 
inco-m-t  isobar. ic  turbulent  boundary  layer  profile 
M  fyi  subjected  to  s~ill  transonic  perturbations 
due  to  an  impinging  weak  normal  shock.  In  the 
prictic-i  Reynolds  number  ran"  of  interest  here 
(Re^  1',’  to  lr)a]  wo  purposely  employ  a  non- 
asymptotic  trible-deck  flow  model  2  m  the  turbu¬ 
lent  boundary  layer  patterned  after  the  Lig^thill- 
Strat ford-Hond >  approach  that  has  proven  highly 
successful  in  treating  a  variety  of  other  problems 
involving  turbulent  boundary  layer  response  to 
strong  r-a.^id  adverse  pressure  gradients  .and  which 
is  supported  by  a  large  body  of  transonic  and 
supersonic  interaction  data.  The  resulting  flow 
model,  rig .  2,  consists  of  an  inviscid  boundary 
value  problem  surrounding  a  shock  discontinuity 
and  underlaid  by  a  thin  shear  stress-d isturoance 
sublayer  that  contains  the  upstream  influence  and 
skin  friction  perturbations.  An  approximate 
analytic  solution  is  further  achieved  by  assuming 
small  linearized  disturbances  ahead  of  and  behind 
the  nonlinear  shock  jump  plus  neglect  of  the 
detailed  shock  structure  within  the  boundary  layer, 
which  give  accurate  predictions  for  all  the  pro¬ 
perties  of  engineering  interest  when  M  -»  1.00.  The 
resulting  equations  can  be  solved  by  operational 
methods  yielding  the  interactive  pressure  rise, 
displacement  thickness  growth  and  the  skin  friction 
behavior  upstream  and  downstream  of  the  shock 
foot.  This  solution  contains  all  the  essential 
global  features  of  the  mixed  transonic  viscous 
interaction  flow  and  detailed  comparisons  with 
experiment  and  Navier -Stokes  numerical 

solutions  have  shown  that  it  gives  a  very  good 
account  of  all  the  important  engineering  features 
of  non-separating  interactions  over  a  wide  range 
of  Mach-Reynolds  number  conditions. 

An  important  and  unique  feature  of  this  inter¬ 
action  theory  is  that  it  employs  for  the  incoming 
turbulent  boundary  layer  velocity  profile  a  very 
general  Composite  Law  of  the  Wall-Law  of  the  Wake 


2.2)  .■.all  Curvature  and  Shock  Obliquity  effects 


Since  SBLI  with  tangential  injection  of ten 
arises  in  flows  on  curved  surfaces,  it  is  desire- 
able  to  account  for  wall  curvature  effects  in  the 
foregoing  interaction  theory.  For  the  small  to 
moderate  curvatures  usually  encountered  (Ki  "  .-2,, 
details  analysis^of  the  transonic  small  disturbance 
flow  in  the  outer  deck  shows  that  while  the 
explicit  new  curvature  terms  in  the  perturbation 
equations  are  of  the  negligible  order  K t  ,  the 
interactive  viscous  displacement  effect  fVom  the 
underlying  decks  eliminates  the  well-known"^ 
inviscid  shoe <  singularity  while  slightly  alt-ring 
the  shock  into  jn  oblique  configuration.  Detailed 
examination  of  the  m  iddle-d'-ck  region  shows  that 
any  new  terms  in  the  inviscid  rotational  di  ;tjr- 
bance  equations  are  of  the  negligible  order  Y-  : 
only  the  curvature  effect  on  the  undi  st  _.i  o'd 
boundary-layer  velocity  and  eddy  viscosity  profile, 
are  of  possible  significance.  Here  again,  the 
explicit  K i  terms  in  the  governing  equations  ot 
this  incoming  flow  j*---  all  negligible;  however, 
curvature  can  moderately  influence  f 10  —  20  *  #  the 
eddy  viscosity  terms,  with  a  consequent,  effect  on 
the  boundary-layer  profile  in  the  form  of  a  skin 
friction  reduction  ind  shape  factor  increase-  j:. 
described  approximately  by  the  relationships® 


a  -  luKf  )  (c 


f  'flat 


a  > 


Li 


■"‘ii'n-t 


C/.i 


where  to  this  order  of  accuracy  the  correspondin'} 
e f i ec t  on  Cf  is  ne  }  1 i g i b 1 y  sma 1 1 .  Note ,  f  or 
example,  that  *  he  typical  value  K6  0.01  yield*; 

a  reduction  md  increase  in  Cf  an§  H/,  of  10  and 

o  1 

0  * ,  respectively.  i'i.e  use  of  Kgs.  (IT  and  (2) 
with  the  Walz  velocity  profile  model  and  K1*.  as  an 
additional  input  parameter  provides  a  good  engin¬ 
eering  account  of  the  moderate  curvature  effects 
on  the  middle-deck  interaction  solution.  Within 
the  very  thin  i nner  disturbance  shear  stress  deck 
it  is  found  yet  again  that  the  explicit  curvature 
effects  on  the  various  inertia,  pressure  gradient, 
and  laminae  viscous  terms  in  the  distur fiance  flow 
equations  are  altogether  negligible.  Moreover, 
because  of  the  extreme  inner -deck  thinness,  the 
eddy  viscosity  curvature  ef f ect^  therei n  can  also  be 
safely  neglected  for  the  high  Reynolds  number  con¬ 
ditions  typifying  most  practical  external  aero¬ 
dynamic  flows. 
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?  C  iCt  ICO  . 


Predictive  results  for  the  typical  influence 
o f  KS  on  SBLI  properties,  which  a  free  with  exper¬ 
imental  observations,  may  :o  found  m  Hoi.  r* ; 
they  show  that  the  curvature  effect  slijntly 
spreads  out  the  interaction,  we.isep.inj  the  adverse 
pressure  gradient  along  th»s  wall,  du*-1  :r  in-..»r  1  ly  to 
the  increased  shape  factor.  Since  the  cat”  .tare 
effect  slight! y  reduces  the  i neon.  *  n- :  no u n- i  : r y - 1» ye r 
velocity  profile  fullness  and  spreads  out  the 
interact  ion,  it  further  acts  to  tmekon  the  down- 
s  t  r  e  am  b«  i  j  n  i  .r  y  1  a  ye  r  while  s  1  i  a  h  1 1  y  l  nc  r  e  a  s  i  n-  i 
the  loe  il  around  the  shock  t'.vf  owing  to  the 
reduced  interactive  pressure  uri-dient.  As  re  i  arris 


the  --light  a  hoc  \  ccliquity  at  tr.e 


e  \  uj  c  a  u  s  ed  :  y  t  h  e  interact  i  v  e  a  i 
■‘•■'ct ,  derail.!  invest  i-uati 
that  it  corrosp.  ndv  with  .«  :ood 


co-diti'n 
sure  r i n e 


o f  do  i  1  ec  t  i or. . 

is  e  ;ui valent  to  ;  norm 


o f  f  ec t i vo  lewe r 


no  and ary  lay.*- r 


:  r-  'Xi;r  ,t  ion  to 

Hence  the  :  r*-.'- 


,  s  i  n  1 9  f 


thereby  allowing  the  onliguity  effect  to  ho  vear- 
»tol'.  .»cco:r..  via  to  i  in  the  r resent  the  ry^  :  f:  : .  A'. 


Extension  to  Include  Tmjentiil  In  tection 


We  recall  that  the  .» forewent  loned  interaction 
theory  feeds  j;an  ••  known  shock  strength  »nd  in 
Assumed  Incot ing  turbulent  velocity  profile  r.odel 
character  iced  by  the  overall  :j  riveter  :}.  ,  Ci  o 

1 1 

and  Re  . + .  In  the  problem  at  hand  we  h  ivo  a  now 
uniqueshape  of  velocity  profile  that  exists  due  to 
tangential  blowing.  In  this  section  we  will  ne 
concerned  wich  modeling  such  a  pro:  i It?  and  its 
associated  wall-region  eddy  viscosity  behavior  ny 
a  convenient  set  of  parameters  tnat  c:i.»r  ictoi  1  ze 
the  essential  new  physical  fo.itjros  and  yet  are 
flexible  er.ouuh  to  sccomod ate  liter  r pec i f i c  data 
ind  to  allow  parametric  sensitivity  studies. 

‘.wien  air  is  tangent i ally  injected  thr-vign  a 
slot  oi  height  h  into  an  overlying  ojundiry  layer 
it  forms  a  jot  which  is  entrained  by  the  surround¬ 
ing  :  low  -Fig.  a  a )  .  Immediately  .lownr.tn?  r  of  the 
si i it,  strong  Tixmj  of  these  flow:  occurs  in  a 
com  pi  lc  ited  sinner  which  may  not  be  validly  treated 
by  t sand  try  layer  theory;  in  any  event,  the  result¬ 
ing  compos i te  velocity  profile  assumes  :  uni. rue 
character  with  a  maximum  and  i  minimum  vl'iu.  Ab>. 

As  thn  flow  proceeds  further  downstream,  experi¬ 
mental  studies  have  shown  that  the  minimum  is 
rapidly  eliminated  by  farther  mixing  so  that  when 
x  '>  h,  tiie  profile  attains  »  f all > -developed 
"jet-bulged"  shape  (Fig.  Ac '  composed  of  an  unblown- 
type  of  turi  ulent  boundary  layer  profile  plus  a 
will  jet  component  c*nt  lining  •  »  velocity  maximum 
n*1  ir  the  wall.  A::,  this  t  ul  ly -dove  loped  shape  con* 
veots  downstrei-,  further  mixing  gradually  de¬ 
creases  ami  spreads  out  the  jet  maximum  (Fig.  ad ) 
until  t lie  bound  ay  layer  ultimately  tends  toward 
in  •  »rd in «ry  monotone  profile  shape  in  which  the 
K’t  component  has  been  completely  eliminated  by 
entrainment.  In  the  present  study  of  weak  tran¬ 
sonic  nor  -  1 1  shock  interaction  with  the  boundary 
1  v/er  downstre tm  of  a  tangential  injection  slot, 
we  will  deal  with  the  case  where  the  shock  inter¬ 
acts  with  a  jot-bulgen  type  of  profile  (Fig.  Ac); 
this  i  •  the  most  interesting  encountered  in 


While  possess  in  s  ■■  boundary  layer  profile 
shape  that  c  in  be  analytically  modeled  in  a  manner 
ip:  rapr i ate  to  the  SBLI  eolation  (see  below),  this 
case  also  permits  t  simplified  treatment  of  the 
eddy  viscosity  aspects  of  the  interactive-  decks  in 
the  boundary  layer,  as  follows.  Experimental 
studios**’*^  have  shown  that  the  usual  Law  of  the 
Wall  sehavior  and  its  associated  mixing  length 
eddy  viscosity  model  applies  to  the  lower  portion 
relow  the  jet  maximum  when  the  injection  effect  is 
3.m  ill-to-moder  *te  (iu^  u  >  l.C).  Since  the 
thin  distur:  jnce-snear  stress  inner  deck  of  the 
SBLI  region  lies  well  within  this  Law  of  the  Wall 
region,  while  there  are  no  eddy  viscosity-associ¬ 
ated  perturbation  term.s  in  the  overlying  middle 
dec*  owing  to  the  inviseid  frozen-turbulence 
nature  of  its  disturbance  flow,  height,  it  can  be 
shown  that  the  form  of  all  the  basic  triple-deck 
Equations  in  the  aforementioned  SBLI  theory  can  be 
carried  over  to  the  present  problem  provided  that 
one  fully  ‘iC counts  for  the  wall  jet  effects  on  the 

•undisturbed  flow  skin  friction  C.  ,  displacement 

*  o 

t::  vs  ness  Roy  no  Ids  number  Re  .  #  and  ( espec  lull  y 
inc.jm.press m-lc  shape  factor  t\:  as  well  as  the  pro¬ 
fib-  distribution  itself.  ^ 

An  appropriate  analytical  model  of  the  incofr- 
m  j  pound  »ry  l  ».er  profile  was  developed  which 
iccojnts  :*.»r  the  essential  new  wall-jet  features 
ol  the  flow  while  also  being  well-suited  to  the 
Li  Jhthill  pressure  disturbance  equation  that  is 
involved  m  the  middle  deck  solution.  It  is 
constructed  a  si  the  sum  of  a  wall- jet  component  and 
«n  "unblown”  component,  where  to  be  consistent 
with  earlier  work  the  latter  is  represented  by 
Wile's  aw  '.*£  the  Will  Law  of  the  Wake  composite 
profile  cr.ur after i?ed  by  the  three  parameters  Re.#f 

d.  c'.d  H.  see  Appendix  A).  Thus  if  y  denotes 
• o  1 1  ' max 

the  he i  ?ht  of  t.he  maximum  velocitv  u  with  .‘u 

max  max 

denoting  the  cor  respond ing  difference  between  u 
snd  tiio  uni  town  velocity  due  to  the  wall  jet 
effect  vCee  Fij.  r» » ,  the  total  profile  is  expressed 


up.-  u  .  (y !  +  .‘.u  (y )  (4) 

WU  1  2 

where  the  will  jet  component  hu(y)  varies  from 
.’era  at  y  =  0  (no  slip)  to  its  maximum  value 
at  y  =  y_  and  then  decays  outwardly 
towards  zero,  becoming  negligible  beyond  some 

characteristic  jet-spread  height  f  .  above  y 

..mix  'max 

iwe  presume  s  .  ♦  v  Above  y  ,  we 

mix  'max  max 

have  followed  tbe  exper imentally-based  work  of 

Carrj,erre  et  .il^and  represented  au  by  a  modified 

sech”  fum:tion  whose  slope  at  y  equals 

- (dp  .  dvJ  such  that  the  ?ofal  composite 

w  a  1 7.  y  - 

profile  correct?-,-  has  a  maximum  at  y  : 

max 


AU 

Au  m  ix 


V  -  >' 


?FCH"  (  (- 


-)  ♦  -M 


SECII  4 


(5) 


•The  regions  upstream  of  the  slot  and  very  far 
downstream  where  the  profile  maximum  has  dis¬ 
appeared  can  ot  course  be  handled  by  the  existing 
"unblown"  vers. ion  of  the  present  SBLI  theory. 


where 


C  C 

»  ,  S [ 1 n ( 1  *  ~)  -  lnll  -  -|)|  (SB) 

is  j  phasing  factor  insurinq  the  maximum  in  total 

velocity  »t  y  and 
max 


C 

x 


(5  .  /AU 
mix  max 


)  I 


«v/l 

max 


is  a  lateral  spreading  constant  (typically  -  .la 

to  avoid  secondary  profile  maxima  above  y  ). 

1  r  max 

Below  y  ,  on  the  other  hand,  we  require  a  func- 
max 

tional  representation  that  gives  a  reasonable 
monotonic  shape  and  matches  smoothly  to  Eu.  (S)  at 
y  .  Furthermore,  we  desire  some  control  over 
tiae  wail  slope  in  order  to  represent  injection 
effects  on  the  local  skin  friction  ACf  .  The 
specific  constraints  on  this  functional  choice  are 
(a)  only  one  maximum  in  the  total  composite  profile 

at  y  *  y  ,  (b)  a  match  with  the  value  and  slope- 
max 

of  the  upper  Au(y)  function  at  y  ,  and  (c)  posi¬ 
tive  values  of  the  non-dimens ionaf  slope 


Id  (Au.'ue) 

d  iy .  y  i 

max  J  ^ 

le  id inq  to  physically  reasonable  skin  friction 
increments 


r 

f 


(6) 


Now  condition  (a)  so  severely  restricts  the  class 
of  monotone  functions  it  admits  that  no  general 
solution  can  be  generated  to  accommodate  a  com¬ 
pletely  arbitrary  combination  of  conditions  (b) 
and  (c);  whit  can  be  found,  however,  are  functions 
which  allow  either  an  arbitrary  choice  of  all 

three  parameters  S  ,  AU  ,  y  within  a  restric- 
w  max  max 

tive  range  or_  the  choice  of  a  wide  range  of  values 
for  the  two  key  parameters  AU  ,  y  with  S 
then  consequently  determined  7&ut  s¥?3fl  within  an 
interesting  range  of  resulting  values) .  One  such 
function  which  has  proven  quite  satisfactory  for 
the  purposes  of  this  investigation  is 


Au 

ue 


c3  |exp(C2y 


1) 


(y  <  y_„)  <7A> 

max 

where  the  aforementioned  matching  conditions  are 
fulfilled  if  the  constants^  2  ^  satisfy  the  three 
simultaneous  relations 


(exp  C2 


1) 


C.  -  C,C  exp  C  -  -  y 
1  3  2  2  max 


C  -  C 
1  3 


(7B) 

(7C) 

(7D) 


This  trio  is  readily  solved  numerically  during  the 
implementation  of  the  velocity  profile  model  by 
using  a  standard  non-linear  simultaneous  root- 
finder  subroutine. 


The  aforementioned  provides  a  smooth, 
piecewise-cont inuous  and  physically  realistic 
analytical  model  of  a  f ully-turbulent  boundary 


layer  downs t re am  of  i  tangential  injection  slot; 
it  captures  the  velocity  over shoot  and  negative 
vorticity  region  features  unique  to  this  kind  of 
flow  (Fig.  h )  while  containing  sufficient  basic 
parameter  i /.»t  ion  to  permit  sensitivity  studies  of 
how  the  jet-bulge  effect  influences  the  SBLJ  zone. 
Moreover,  it  has  the  advantage  of  allowing  current 
and  later  experimental  data  on  turbulent  wall-jet 
boundary  layer  behavior  to  be  incorporated  into  th 
interaction  study  without  tying  the  present 
research  down  to  the  much  more  difficult  and 
lengthy  effort  of  such  experimental  studies.  The 
weak  boundary  layer  compressibility  effects  on 
this  profile  tor  adiabatic  transonic  flow  are 
quite  satisfactorily  accounted  for  by  the  referenc 
temperature  method.1® 


3.2)  Implementation  of  the  Intended  Theory 


The  foregoing  approach  may  be  implemented  by 
several  straightforward  modifications  to  the 
existing  computer  program  for  the  zero-blowing 
SBLI  theory,  as  follows.  To  include  smull-to- 
moderate  wall  curvature  effects  (K6  <  .01),  wo 

add  Ki  as  an  independent  input  variable  and 
accordingly  modify  the  input  values  of  lb  and  CfQ 

according  to  Eqs.  (1)  and  (2);  furthermore,  we 
eliminate  the  inviscid  curvature  singularity, 
altering  the  normal  shock  to  a  slightly  oblique 
one  at  the  boundary  layer  edge,  by  modifying  the 
input  effective  normal  shock  Mach  number  according 
to  Ft).  (3).  The  influence  cf  tangential  injection 
is  accomodated  by  introducing  the  two  new  input 

parameters  AU  U  and  y  /6  ,  characterizing 
max  e  7  max  o 

the  magnitude  and  height,  respectively,  of  the 
wall  jet  component  effect;  in  addition,  values  of 

the  auxiliary  parameters  C  and  S  can  be  set 
4  x  w 

within  certain  restricted  ranges.  The  program 
subroutine  which  evaluates  the  Walz  turbulent 
boundary  layer  velocity  profile  model  is  modified 
to  add  the  matched  upper  and  lower  wall  jet-com¬ 
ponent  increments  pertaining  to  these  inputs 
(Eqs.  4-7),  using  u  Hefercnce  Temperature-Method 
compressibility  correction  of  the  appropriate 
parameters.  Figure  7  illustrates  some  typical 
boundary  layer  velocity  profiles  containing  these 
tangential  injection  effects.  Using  the  adiabatic 
temperature-velocity  relat ionship 


2  2 

T  =  T  +  (T  -  T  )U  /lb 

W,  AD  e  W,AD  <? 


(8) 


the  associated  Mach  number  profile  M  (y)  and  its 
derivative  dM  /Jy  (which  are  both  needed  in  the 
subsequent  SB?I  solution  routine)  are  calculated, 
the  corresponding  mass  flow  and  momentum  defect 

distributions  1  -  — and  (1  -  are  also 

oeUe  peUe  ue 

integrated  across  the  boundary  layer  to  obtain  the 
values  of  5*/<$  and  0V6,  respectively,  associated 
with  the  wall  jet  effect.  The  resulting  values  of 
the  displacement  thickness  and  shape  factor  are 
shown  in  Figs.  8A  and  8B,  to  illustrate  how  the 
mass  and  momentum  addition  to  Lhe  boundary  layer 
from  the  wall  jet  substantially  decreases  6*  and 
produces  a  greater  profile  **£ ullness"  +  ref lected  in 
a  siqnif icantly  reduced  shape  factor.  Increasing 

*  It  is  formally  possible  to  obtain  negative 
6*  and  0*  for  sufficiently  large  injection  rates 
<AUmaX»  say);  consistent  with  our  other  assumptions 
however,  we  exclude  such  cases  from  this  study. 


* 


the  height  of  the  jet  maximum  is  seen  to  have  a 
similar  effect,  because  this  enhances  the  effective 
strength  of  the  injection  effect  on  the  boundary 
layer  profile.  Awareness  of  these  overall  integral 
property  effects  proves  helpful  in  interpreting  the 
predicted  interaction  properties  given  below. 

Implementation  of  these  wall  jet-modif ications 
is  quite  straightforward,  except  to  note  that  feed¬ 
back  of  the  aforementioned  modified  integral  pro¬ 
perties  into  the  solution  seguence  must  be  properly 
phased:  since  the  wall  jet  effect  on  the  incoming 
boundary  layer  profile  shape  is  already  included  in 
the  M  (y)  distribution  used  in  solving  the  Light- 
hill  interactive  pressure  equation,  the  feedback 
must  be  done  after  this  pressure  disturbance 
solution  is  carried  out.  Subseguent  use  of  the 
jet-altered  values  of  5*  and  CfG  then  further 
influences  the  local  interactive  displacement 
thickening  and  skin  friction  solution  results.  To 
illustrate  the  importance  of  this  proper  feedback 
of  the  jet-influenced  profile  integral  properties 
a  typical  set  of  interactive  pressure,  displacement 
thickness  and  skin  friction  distributions  predicted 
by  the  aforeme  tioned  extended  theory  are  presented 
in  Fig.  9,  showing  the  various  relative  effects  of 
tangential  injection  compared  to  the  zero  blowing 
case.  It  is  seen  that  the  increased  boundary 
layer  profile  fullness  and  shape  factor  reduction 
due  to  injection  causes  a  significant  str^miwise 
contraction  of  the  interactive  pressure  rise;  this 
is  in  agreement  with  experimental  obsemttonr, 

(see,  e.a ,  Fig.  116,  p.  1123  ol  Ref. 17  ).  Accom¬ 
panying  this  contraction  of  the  interaction  zone, 
the  two  main  effects  of  injection  on  the  ratio 

are  seen  to  act  with  opposite  and  nearly 
egual  influence:  while  the  profile  shape-factor 
effect  of  injection  reduces  A**,  the  cor  responding 
reduction  of  50*  is  approximately  of  the  same  mag¬ 
nitude  so  that  the  overall  change  in  a5V60*  is 
small.  This  implies  that  the  net  injection  effect 
on  .15*  scales  approximately  with  the  correspond¬ 
ing  effect  on  o  *.  Turning  to  the  interactive 
skin  friction  behavior  typified  in  Fig.  9c,  it  can 
be  seen  that  the  increased  Cf  level  due  to  the 
wall  jet  effect  dominates  most  of  the  interaction 
zone  both  fore  and  aft  of  the  shock  except  in  the 
vicinity  of  the  shock  foot;  in  this  foot,  region, 
the  Cf  reduction  due  to  the  steepened  interactive 
pressure  gradient  caused  by  injection  becomes  the 
dominant  effect  and  the  local  value  of  Cf  .  is 
actually  reduced.  Stated  another  way,  the  SBLI 
effect  adversely  counteracts  the  otherwise  favor¬ 
able  Cf  increase  due  to  injection. 

The  »forement ioned  tangential  injection 
effects  on  SBLI  may  be  readily  understood  from  the 
overall  shape  factor  ind  displacement  thickness 
effects  shown  in  Fig.  8:  the  reduced  H  and  6* 
imply  a  thinner  incoming  turbulent  boundary  layer 
with  a  somewhat  higher  Mach  number  deep  in  the 
layer  and  a  fuller  profile  shape  typical  of  a 
favorable  ups tr earn  pressure  gradient  history,  which 
in  view  of  the  demnnstr  ited  sensitivity  of  SBLI  to 
the  shape  factor  H-’ig.  I)  have  the  effect  of 
reducing  the  streimwise  scale  and  interactive 
thickening  while  increasing  the  corresponding  local 
pressure  gradient. 


3.3)  Imbedded  Regions  of  Negative  Vort ic ity  and 
Supersonic  Flow  in  the  Boundary  Layer 

It  has  been  seen  that  the  wall  jet  effect 


results  in  a  strata  of  negative  vorticity  flow 
above  the  maximum  deep  down  in  the  incoming  boun¬ 
dary  layer  profile  (Fig.  6).  Now,  some  earlier 
basic  studies  of  shock  interaction  with  idealized 
shear  flows  isimple  velocity  discontinuities)  sug¬ 
gest  that  such  a  strata  of  vorticity  sign  reversal 
might  significantly  alter  the  character  of  the 
shock  transmission  and  reflection  across  it,  in 
turn  implying  possible  difficulty  with  the  numeri¬ 
cal  solution  across  this  strata  of  the  Lighthill*® 
interactive  pressure  disturbance  equation  in  the 
present  SBLI  theory  (which  involves  a  term 
-  -*P/dV  *  dU0/dy).  We  therefore  examined  this 
point  carefully,  with  the  following  reassuring 
conclusion:  provided  that  reasonable  care  is  taken 
to  insure  high  numerical  accuracy  with  an  appro¬ 
priately  smaller  step  size  iv ,  the  Lighthill  equa¬ 
tion  solution  is  guite  regular  for  any  smooth 
albeit  rapid  variation  in  sign  (dMo/dy)  across  the 
strata.  Hence  the  overall  interaction  solution  is 
modified,  but  not  fundamentally  altered,  by  the 
presence  of  the  negative  vorticity  due  to  the  wall 
jet  effect  and  this  is  straightforwardly  accounted 
for  by  our  modified  velocity  profile  model  in  the 
Lighthill  equation  and  by  the  associated  change  of 
the  integral  parameters.  The  underlyin']  reason 
for  this  lack  of  difficulty  with  rapid  local  vari¬ 
ations  in  either  magnitude  or  sign  of  dMQ/dy  may 
be  found  from  an  analysis  of  the  large  scale 
features  of  Lighthill’s  equation,  which  reveals 
that  its  solution  essentially  depends  only  on 
integrals ,  rather  than  on  local  details,  of  the 
Mo(y)  distribution  across  the  boundary  layer. 

The  presence  of  a  local  velocity  maximum  deep 
within  the  l>oundary  layer  also  raises  another 
possible  dil 1 iculty ,  when  the  wall  jet  effect  is 
sufficiently  large,  associated  with  the  existence 
of  a  strata  of  locally  supersonic  flow  astride  the 
velocity  maximum  (Fig.  10).  When  this  occurs,  it 
is  seen  that  there  are  two  special  cases  where 
dM0/dy  vanis-.hes  at  a  sonic  point  within  the  boun¬ 
dary  layer  and  where  a  local  transonic  singularity 
in  the  Lighthill  pressure  equation  solution  there¬ 
fore  will  occur:  (a)  at  a  tangential  injection 
rate  where  Umax  just  goes  sonic,  and  (b)  at  a 
slightly  higher  rate  where  the  local  minimum  U  goes 
sonic  higher  up  in  the  boundary  layer.  In  these 
two  isolated  cases,  there  is  a  local  breakdown  of 
the  linear i zation  underlying  the  Lighthill  equation 
and  the  resulting  transonic  singularity 
which  causes  fundamental  difficulties  with  the 
numerical  solution  of  this  equation  that  can  only 
be  cured  by  restoring  (at  least  locally)  the  appro¬ 
priate  non-linear  transonic  correction  term.  For 
all  other  maximum  wall  jet  velocities  (including, 
interestingly  enough,  the  so-called  "overblown" 
cases  where  Mmax  >  Mg),  the  boundary  layer  con¬ 
tains  only  one  local  sonic  point  that  is  well- 
removed  from  dM^^/dy  =  0  (for  subsonic 

Umax  it  lies  above  ymax  while  for  supersonic  Umax 
it  lies  below).  in  such  normal  cases,  no  funda¬ 
mental  difficulties  were  discerned. 


4.  Discussion  of  Parametric  Study  Results 

The  present  theory  has  been  used  to  carry 
out  a  systematic  study  of  how  the  key  tangential 
injection  parameters  influence  the  essential  pro¬ 
perties  of  a  subsequent  SBLI  zone.  We  now  present 
and  discuss  the  results. 


5 


4.1 


Interactive  Pressure  and  Displacement:  Thick¬ 
ening 


Typical  pressure  distributions,  showing  the 
strong  systematic  contraction  of  the  streamwise 
interactive  scale  with  increasing  strength  of  ♦‘he 
wall  jet  component -ef feet ,  are  illustrated  in 
Figure  11.  A  comprehensive  summary  of  such 
results  showing  the  upstream  and  downstream  influ¬ 
ence  distances  (the  distance  ahead  and  behind  the 
shock  where  the  pressure  rise  is  5*  and  95*, 
respect i vely ,  of  the  overall  shock  jump  value)  ire 
presented  in  Figures  12  and  13  as  a  function  of 
both  the  magnitude  and  location  of  the  jet  velocity 
maximum  for  a  typical  supercritical  flow  of  M  - 
1.20,  Additional  plots  showing  the  influence  of 
the  incoming  (unblown)  shape  factor,  shock 
strength  and  Reynolds  number  on  the  wall  j*»t 
effects  are  presented  in  Figures  14-17.  Taxon 
overall,  these  results  show  that  tangenti  il 
injection  can  significantly  reduce  the  over  *11 
upstream  influence,  and  strongly  reduce  the  down¬ 
stream  streamwise  scale  of  the  interaction  to  .> 
degree  comparable  to,  or  greater  than,  the  unblown 
shape  factor  and/or  Mach  number  effects.  When 
non-di"ensionalized  in  terms  of  ^  ,  the  results 
are  not  very  sensitive  to  Reynolds  number . 


predicted  slightly  'Fig.  20a):  when  the  shock 
obliguity  effect  is  modeled  as  described  above,  it 
is  seen  that  th"  present  theory  and  experiment  are 
in  jo'xl  t  )r**or  r-r.*  over  a  wide  range  of  Reynolds 
n  ur  i  >e  r  s  . 

Turning  to  the  effect  of  injection,  we  first 
nut-*  l  r or  the  typical  behavior  of  the  inter  active 
Cf  diatr  1/ jt tor.  around  the  shock  foot  (see,  e.g,, 
rig.  9c,  that  the  net  effect  is  expected  to  de- 
ere  me  Cf  n  (notwithstanding  the  overall  upstream 
and  downstream  increase  in  Cf  otherwise  due  to 
injection;.  As  shown  in  Fig.  21,  this  is  indeed 
tojnd  t-,  be  t!v-  case:  the  wall  jet  effect  of 
increasing  the  1-ical  interactive  pressure  gradient 
is  seen  to  hasten  the  onset  of  incipient  separation 
it  the  shock  foot  for  a  given  Reynolds  number  flow; 
in  the  sens**  that  separation  occurs  at  a  slightly 
lower  shock  number  as  is  increased.  This 

i*  of  course  in  sharp  contrast  to  the  well-known1^ 
favorable  effect  of  injection  in  delaying  separa¬ 
tion  observed  for  purely  subsonic*!  flows  with  a 
prescr ii»f*d  adverse  pressure  gradient,  and  is  due  to 
the  fact  that  the  interactive  pressure  gradient 
enhancement  effect  of  tangential  injection  in 
locally  reducing  Cf  is  absent  in  the  latter  flows. 


The  corresponding  systematic  influence  of 
injection  on  the  rclat i ve  interactive  displacement 
thickness  distribution  i5*(x)/40*  is  illustrated 
in  Figure  Id,  where  we  see  that  the  effect  on 
iF'(x)  and  50*  largely  cancel  over  a  wide  range  of 
wall  jet  strengths  when  presented  in  this  ratiood 
manner.  However ,  there  is  a  significant  injection 
effect  on  the  streamwise  slope  of  a*(x)  at  tfm 
shock  foot,  which  relates  to  the  effective  "vis¬ 
cous  wedge"  angle  sensed  by  the  outer  inviscid 
flow;  this  effect  is  illustrated  in  Fig.  19,  where 
the  strong  increase  of  this  slope  with  wall  jet 
strength  may  he  clearly  seen. 


4.3;  Down st re  «r  Effects 
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4.2)  Incipient  Separation 

The  present  theory,  although  it  breaks  down  at 
separation,  does  yield  a  useful  indication  of 
incipient  separation  where  Cf  -*0,  owing  to  the 
particular  attention  paid  to  Tfte  treatment  of  the 
local  interactive  skin  friction  behavior.  Since 
this  indication  is  of  great  practical  interest,  a 
parametric  study  of  incipient  separation  conditions 
inherent  in  th**  present  theory  was  carried  out. 

As  a  basis  for  comparison,  the  results  for 
flow  without  any  tangential  injection  are  shown  in 
Fig.  20a  where  the  shock  Mach  number  above  which 
incipient  separation  occurs  is  plotted  as  a  func¬ 
tion  of  the  Reynolds  number  with  the  shape  factor 
as  a  parameter;  also  shown  in  the  figure  is  the 
approximate  experimental  boundary  determined  by  a 
careful  examination  of  a  large  number  of 

transonic  interaction  tests,  besides  Nussdorfer's 
M  -  1.30  criterion  for  turbulent  flow.  It  is 
seen  that  the  theoretical  prediction  of  a  gradual 
inert  ju.  in  the  incipient  separation  Mach  number 
value  with  Revnolds  number  is  in  agreement  with 
the  trend  of  this  data.  The  theoretical  predic¬ 
tion  of  only  a  small  influence  of  shape  factor  on 
the  incipient  separation  conditions  is  also  borne 
out  by  more  recent  dat as  indicated  in  Fig.  20b. 

We  note  here  that  the  absolute  values  of  the 
incipient  separation  Mach  number  predicted  on  the 
basis  of  a  normal  shock  are  consistently  under¬ 


flow,  wo  have  soon  above  that  tangential 
injection  reduces  the  SBLl  displacement  thickness 
growth  while  increasing  the  downstream  post-inter¬ 
active  Cf  ( Fig .  9c).  Alternatively,  we  may  view 
SBLI  as  increasing  the  downstream  rf * ,  and  hence 
counteracting  the  thinning  effect  otherwise 
obtained  by  the  wall  jet,  while  reducing  the 
injection-produced  Cf  enhancement;  both  these  S13I.J 
effects  make  the  boundary  layer  less  resistant  to 
separation  in  any  subsequent  adverse  pressure 
gradient  region  it  may  encounter,  md  hence  dimin¬ 
ish  the  effectiveness  of  injection  in  otherwise 
delaying  downstream  separation.  Regarding  the 
skin  friction,  these  conclusions  are  summarized  in 
Fig.  23,  where  there  is  shown  the  typical  influ¬ 
ence  of  increasing  wall  jet  strength  on  the  |>or.t- 
inter active  Cf :  it  is  seen  that  while  weak  injec¬ 
tion  at  first  increases  it  slightly  due  to  the 
corresponding  increase  in  CfQ,  stronger  injection 
rates  have  the  opposite  effect  of  lowering  it  (as 
well  as  Cfmin)  because  of  the  intensified  adverse 
pressure  gradient  effect. 


5.  Concluding  Remarks 

Viewed  overall,  the  present  study  has  shown 
that  the  usual  favorable  tangential  injection 
effects  of  thinning  out  and  delaying  the  separa¬ 
tion  of  turbulent  boundary  layers  in  subsonic 
flow  can  be  significantly  compromised  by  transonic 


shoe*,  boundary  layer  interaction.  Conversely,  such 
injection  was  :,een  to  appreciably  reduce  the 
streamwise  extent  of  an  SBL1  zone  ilbei  t  with  the 
allied  consequence  of  intensifying  the  local  inter¬ 
active  adverse  pressure  gradient  and  onset  of 
shoe*  foot  se;'-'*  jtion.  It  has  further  been  estab¬ 
lished  that  a  fund amenta lly-oased  triple-deck 
theory  of  SBLI  with  injection  is  r.ow  available  to 
treat  these  effects  in  either  extern. >1  or  internal 
supercritical  flow  field.:;  moreover,  this  theory 
has  been  constructed  to  serve  as  a  locally  ins^rt- 
m  le  interactive  module  istride  the  mviscid  shock 
location,  driven  by  the  attendant  local  boundary 
layer  properties  including  an  arbitrary  non- 
equ  1 1  i:  r  ium  shape  factor.  Consequently  it  would 
;e  possible  to  investigate  in  the  future  interest- 
mu  problems  of  blowing  in  superer i t ical  flow 
:  ields,  including  the  use  of  tangential  in  jection 
to  modifv  the  influence  of  SBLI  upon  the  viscous 
trail  mu  e-die  ertect  ot  supercritical  airfoils*-**' 
:.nd  the  me  Iasi  n  of  SBLI  effects  in  viscouc- 
mviccid  f  low  r  ield  analysis  pr  oar  airs  for  circulu- 
t  ion-cont  rv-l  led  airfoils  and  wings  flying  .«t 
j uporcr it icil  : liaht  speeds. 
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Appe nd i x 

Because  of  its  convenient  analytical  for:", 
s.vunte  denied  representation  of  the  cord  mod 
Lav;  of  the  Wall  -  Law  of  the  Wake  behavior  and 
toner ality,  we  have  adopted  W.ilz'a  model’  for  the 
incoming  tur:  ulent  sound  ary  l  ayer  np.n  ro.r  of  the 
intersection.  for  the  low  Mach  number  small  heat 
transfer  conditions  rpropriute  to  transonic  inter- 
lotion ,  it-  a  y  :  e  s . » 1 1  •  factor  1 1  y  oo  r  r^ct  ed  i  o  r 
J  •  i  e :» si:  lilt y  e  t  f -vt  n  ivy  t  he  ho ker  t  Re f  e r once 
?e:r  n  ttur^  ••thod  w:\ich  under  these  conditions 
is,  m  *  act ,  oevpar  iblc  m.  accuracy  to,  but  fir 
aim:  ler  to  implement  tl.  in,  the  Van  Priest'**  com¬ 
press  ibi  1 1 1  y  t ions  format  ion  approach. 

Let  -  be  doles'  | incompressible )  Waxo  Function 
y  s  and  denote  for  convenience  R  .41  Re  5* 

III  *■  i  iT.T  ^*‘j  with  ..  -  .76  and  >  ■-  1.4  £or  a 
perfect  gast  tnen  the  compressible  form  of  Walz's 
composite  profile  may  be  written 


lT  ’  '  ' 1  *  H 


1*R-  ,  - 3Rn. 

•  n  (- — — 1  -  [.21)  ♦  .»>  *>R")e  I 


subject  to  the  following  condition  linking  JT  to 
A  ind  Re'* 


21-  *  .nil*Ri 


*  T  . 

~r  {^r) 

2  T 


Eqs .  (A-l)  and  (A- 2)  have  the  following  lies ir able 
properties:  (a)  for  n  N  . 1 1 j  or  so  Uc  Up  is  domin¬ 


ated  by  a  Law  of  the  Wake  behavior  which  correctly 
satisfied  both  the  outer  limit  conditions  L1  V  +  1 
and  dU  ,  dy  ■»  0  as  n  -*  1 ;  (b)  on  the  other  hand,  for 

ver v  small  values,  V  assumes  a  Law  of  the  Wall  - 

o 

type  behavior  consisting  of  a  logarithmic  term  that 
is  exponentially  damped  out  extremely  close  to  the 
wall  into  a  linear  laminar  sublayer  profile 
U  U  -  Rn  as  ♦  0;  (c )  Eq.  (A-l;  may  be  differen¬ 
tiated  w.r.t.  n  to  yield  an  analytical  express:  n 
for  d’J  dy  also,  which  proves  advantageous  in  solv¬ 
ing  the  middle  and  inner  deck  interaction  problem::, 
(see  text i  where  dM  .  dy  must  be  known  ard  vanish  at 
the  boundary  layer  edge. 

The  use  of  the  incompressible  form  of  Lq .  (A-l; 
in  the  defining  integral  relations  for  5  *  and  b  * 
yields  the  following  relationship  that  links  the 
w.jke  parameter  to  the  resulting  incompressible  shape 


(l  +  1 .  ~  9 " 
1  ♦ 


f'qs.  ( A—  2 j  and  (A- 3)  together  with  the  defining 
relation  for  R  enable  a  rather  general  and  conven¬ 
ient  parameter izat ion  of  the  profile  {and  hence  the 
interaction  that  depends  on  it1  in  terms  of  three 
important  physical  quantities:  the  shock  strength 
(Me.),  the  displacement  thickness  Reynolds  number 
Re5*,  the  wall  temperature  ratio  T^T  and  tne 
shape  factor  Hj  that  reflects  the'prior  upstream 
history  of  the  incoming  boundary  layer  including 
possible  pressure  gradient  and  surface  mass  trans¬ 
fer  effects.  With  these  parameters  prescribed,  the 
aforementioned  three  equations  may  be  solved  simul¬ 
taneously  for  the  attendant  skin  friction  C7  ,  the 

1  o 

value  of  R  and,  i i  desired,  the  -  value  appropriate 
to  these  flow  conditions. 
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